Excited states of a dilute Bose-Einstein condensate in a harmonic trap 
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The low-lying hydrodynamic normal modes of a dilute Bose-Einstein gas in an isotropic harmonic 
trap determine the corresponding Bogoliubov amplitudes. In the Thomas-Fermi limit, these modes 
have large low-temperature occupation numbers, and they permit an explicit construction of the 
dynamic structure function ^(q, tj). The total noncondensate number A'''(0) at zero temperature 
increases like R^, where R is the condensate radius measured in units of the oscillator length. The 
lowest dipole modes are constructed explicitly in the Bogoliubov approximation. 
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I. INTRODUCTION 

Recent low-temperature experiments have observed Bose-Einstein condensation in alkali atoms confined to har- 
monic potentials Subsequent investigations have detected the first few low- lying collective modes of the con- 
densate [^||. These results have stimulated a great deal of theoretical activity pertaining to trapped dilute Bose 
systems. Most of this work relies on the Bogoliubov approximation which assumes that only a small fraction of 
the particles are excited out of the condensate. Clearly such a description fails completely at the onset temperature 
To for Bose condensation, where the condensate occupation number N(){T) vanishes, but it does provide a very useful 
description of the low-temperature behavior of a dilute condensed Bose gas. 

In a typical experiment of Refs. [|l],^, the condensate is "large," in the sense that its diameter Rq far exceeds the 
characteristic width do = \JJiJmLjo of the single-particle ground state of the harmonic trap. (Typically, do ~ 1 /xm for 
the traps of Refs. [|l|,|.) This result can be easily seen by comparing the relative contributions of the kinetic energy, 
the trap potential energy, and the interaction potential energy to the total energy of the gas. The mean kinetic 
energy (r)o ~ Nq(Ji' jlml^ diminishes as the condensate grows, whereas both the mean single-particle confinement 
energy (y)o and the mean two-particle interaction energy (i.e., the Hartree energy) (yf/)o increase with condensate 
size. (The subscript denotes an expectation value in the selfconsistent condensate.) For a large condensate, the 
kinetic energy is then much smaller than both the potential energies. Neglecting the kinetic energy in comparison 
with the two potential energies gives rise to the so-called "Thomas-Fermi" (TF) approximation, which provides a 
simple description |^ of the spatially varying condensate density no(r). 

The two-body repulsion between atoms can be characterized by an s-wave scattering length a > 0, or equivalcntly 
an s-wave pseudopotential g = inh^a/m. The characteristic Hartree energy of the condensate is then {Vh)o ^ 
gN^/RQ ~ h^aNQ/mRQ. In the Thomas-Fermi limit, this Hartree energy must dominate the kinetic energy (T)q. 
Their ratio defines an important dimensionless quantity 

Noa {VH)a 
da (T)o 

The TF limit therefore requires not only a macroscopically occupied condensate (iVo S> 1), but also the more stringent 
condition ||] t^q = Noa/do ^ 1. 

In the Thomas- Fermi limit, the characteristic radius Rq of the condensate is determined by the balance of 
{Vh)o ^ h^aN^ /mRQ and the trap potential (V)o, which varies as ^ NQ{muj'^RQ) for a harmonic confining po- 
tential. Minimizing {V -\- Vh) with respect to i?o yields A^o ^ ^o/('^o'^)- Introducing the dimensionless condensate 
radius R = Ro/do, we find 770 ~ R^ in the Thomas-Fermi limit. 

The collective spectrum of a confined cloud of Bose condensed atoms differs qualitatively from that of a droplet of 
liquid helium of comparable dimensions. The helium droplet has nearly uniform density except near a thin surface 
layer. Its bulk collective modes are simply the compressional sound modes of a uniform fluid, whose quantization is 
set by the boundary conditions at the surface of the drop. For a drop of radius Rq, the minimum wavenumber of a 
phonon varies as I/Rq, and the energy spacing between the collective modes is therefore of order Ks/Rq, where s is 
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the bulk speed of sound. As the drop size grows, the minimum excitation energy tends to zero for a large drop, and 
the spectrum approaches the gapless spectrum of bulk helium, as expected. 

In a confined compressible atomic gas, however, the density of the condensate varies with the size of the cloud, in 
contrast with the nearly constant density of a drop of liquid. We have seen that, in the Thomas-Fermi limit, a cloud 
of radius i?o contains A^o ~ R^{do/a) condensed atoms. The mean density of the condensate is then Nq/Rq ^ E? . 

The speed of sound in a dilute Bose gas varies as the square root of its density so that s ^ R. The typical 
energy spacing between quantized sound modes in a harmonically confined cloud then varies as s/ R ^ R^ . Thus while 
the maximum wavelength grows linearly with the size of the cloud, the increased density of a larger cloud raises the 
speed of sound proportionally, and the minimum excitation frequency is independent of the radius of the cloud . 

The excited states of spatially inhomogeneous Bose condensates can be studied by two equivalent approaches. The 
original work of Bogoliubov Q emphasized the underlying quantum character of the problem. His treatment of 
the uniform condensate is readily generalized to the nonuniform case to yield the Bogoliubov equations |9[[Tc[| that 
constitute the quantum-mechanical Schrodinger equations for a pair of coupled amplitude functions. The eigenvalues 
of this problem are the energies of the elementary excitations of the Bose condensate; the corresponding amplitudes 
determine the spatially varying noncondensate density. 

More recently, a hydrodynamic approach has proved valuable in determining both the excitation frequencies and 
the normal- mode amplitudes for the low-lying excited states of a Bose condensate in a harmonic trap in the Thomas- 
Fermi limit [ pT| . It is not difficult to prove that these two descriptions are completely equivalent in the Bogoliubov 
approximation |p^ , p^ . The present work exploits this feature to determine the contribution of the low- lying collective 
modes to the zero-tcmperature noncondensate density. In particular, the occupation number of these low-lying modes 
is large in the TF limit, and a reasonable cutoff for the sum over all such modes suggests that the total noncondensate 
number N' then scales as R^ in this limit. Since scales with i?^, the Bogoliubov approximation that N' <C Nq 
necessarily fails for sufhciently large R (or N). 

We review the basic formalism in Sec. II and summarize the equivalence between the Bogoliubov and hydrodynamic 
descriptions in Sec. Ill, along with the physical properties of the noncondensate. The TF solution for the hydrodynamic 
normal modes is reviewed briefly in Sec. IV. In Sec. V we use the resulting eigenmodes to determine the corresponding 
Bogoliubov amplitudes and low-temperature noncondensate occupation. Finally, in Sec. VI, we use the eigenmodes 
to construct the dynamic structure function S{(\,uj) in the TF limit. An Appendix contains an explicit construction 
of the exact lowest dipole modes of an interacting Bose gas as well as the corresponding modes in the Bogoliubov 
approximation. 



II. BASIC FORMALISM 



We briefly review the Bogoliubov approximation for a nonuniform condensed Bose gas, which was introduced 
independently by Gross |l4| and Pitaevskii to study vortices and their excitations. (In this context, the condensate 
has a uniform density no except in the immediate vicinity of the vortex core.) 

At low densities, the two-particle interaction potential may be replaced by a short-range pseudopotential with 
Viv) w g(5(r), where g is expressed in terms of the s-wave scattering length a through the relation g = 47ra7i^/m [ p^ . 
The present work considers only "repulsive" interactions with 5, a > 0. 

For a uniform condensate with density riQ, a small deformation of the condensate wave function with spatial scale 
A involves a squared gradient {i.e., kinetic) energy ~ ?i^/2mA^. It is useful to define a "coherence" (or "correlation") 
length 



through the balance between this kinetic energy and the repulsive interparticle potential energy gn^. The coherence 
length becomes arbitrarily large for an ideal Bose gas {i.e., when a —> 0). 

In the Bogoliubov approximation, the fractional depletion of the condensate N' /N is of order [^,0 \/nQa^. For 
the Bogoliubov approximation to be valid, N' /N must be small, which implies the number of particles per scattering 
volume uqu'^ is much smaller than unity. (This condition is strongly violated in liquid ^He.) Equivalently, we require 
that ^/a ~ (ripa^)"^/^ be much greater than unity. 

The presence of an external trap introduces another length Rq that characterizes the spatial size of the condensate. 
If ^ ^ R(), then the system resembles an ideal Bose gas with negligible interactions. If ^ ^ Rq, however, the system 
differs qualitatively from an ideal Bose gas. In the experiments on Bose-condensed sodium atoms in Ref. [H, the 
scattering length is a ^ 4.9 nm, the trap has a characteristic oscillator length do ~ 1.9 /im, and there are Aq ^ 5 x 10^ 
condensed atoms. The equivalent mean (isotropic) condensate radius Rq ^ 20 /im implies a central condensate number 
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density ng ~ 4 x 10^° m~^. Then ^ 0.14 /im ^ Rq, and the system is indeed a dilute interacting Bose gas with 
a ^ ^ <^ Rq, definitely far from ideal. 

For simplicity, we consider here only a spherical harmonic trap, with 

F(r) = ^mcuy, (3) 

and a characteristic length scale 



do = d (4) 

V mojQ 

corresponding to the Gaussian width of the single-particle ground state of a single particle of mass m in the trap. When 
the trap contains a large number A^o of condensed particles, their mutual repulsion causes the cloud of atoms to expand. 
The actual condensate density ^ Nq/R^ is then much smaller than the simple estimate ~ iVo/dg. Specifically, when 
the Thomas- Fermi parameter 770 = Noa/do is sufficiently large, the dimensionless radial expansion factor R = Ro/do 
is H] (15?7o)^/^. This reduction in the particle density (by a factor of order 1]^ '^^^) means that the system remains 

dilute for rj^^ [a/ d^Y <C 1. Put another way, A^'o must be much less than (da/a)^, which is ~ 10^^ for Bose- condensed 
sodium atoms in Ref. p^ ]. 

A spatially nonuniform Bose condensate is characterized by a condensate wave function ^'(r) that can be normalized 
to the total number of condensate particles, i.e., J d'^r j^fp = A'o. Then no(r) = |^'(r)P is the condensate particle 
density. For a dilute Bose gas at low temperature, ^! obeys a selfconsistent nonlinear Schrodinger equation known as 
the Gross-Pitaevskii (CP) equation 

{T + V + VH-fi)'^ = 0, (5) 

where T = — 7i^V^/2to is the kinetic-energy operator, V{r) is the trap potential energy operator, Vh(j) — g\"^{r)\'^ — 
gno{r) specifies the mean (Hartree) pseudopotential due to the condensate, and fj, is the chemical potential. For 
a stationary condensate, ^' can be taken as real, but the generalization to a complex condensate wave function 
x^) = e*"^ with superfluid velocity Vj = (?i/m) VS* is not difficult. (Such a complex condensate could describe, for 
example, a vortex [p|jl2[mjl5|| . ) 

If the mean condensate kinetic energy (T)o = / d^r'i'* is negligible compared to the potential energies (Vff)o 
and (V^)o @, which holds for 770 ^ Ij then the Gross-Pitaevskii equation [Eq. (^)] can be approximated by its last 
three terms. In this limit, for each spatial position r, either the condensate wave function vanishes or the condensate 
density satisfies the "Thomas-Fermi" approximation 

VH{r) = gm{v) = [^i - V{v)\ e[^x - V{v)\ , (6) 

where 0(x) denotes the unit positive step function. For a spherical harmonic trap with oscillator length dg and 
oscillator frequency lhq, the TF condensate density is an inverted parabola that vanishes beyond a dimensionless 
cutoff radius R w (15r/g)-'^/^ defined in terms of the chemical potential /i = ^TiloqR^ 

The validity of this TF approximation has been investigated both numerically for various values of the dimensionless 
parameter R ||l9|,^ and analytically through an expansion in powers of the "small" parameter This latter 

treatment shows that the TF approximation fails in a thin surface region of thickness ~ do(dg/i?g)^/'^ ^ dgi?"-'^/'^, 
where the formally negligible correction terms eliminate the singularity of \d^'YY/dr\^ at the condensate surface. As 
a result, the condensate kinetic energy acquires a logarithmic correction of order (fiuJo/2R'^) \nR. 

In an ideal Bose gas at zero temperature, all the particles are condensed in the single-particle ground state of 
the trap. Repulsive interactions excite a (small) fraction N' /N of the particles out of the condensate, even at zero 
temperature. These excited particles occupy the various normal-mode eigenstates that satisfy the (linear) Bogoliubov 
equations |J|,|l^,|l|,[| 

Cuj — VhVj — EjUj, (7a) 

— VnUj + Cvj = —EjVj, (7b) 

for the coupled eigenfunctions Uj{v) and Wj(r), and the associated eigenvalues Ej. Here, j denotes a complete set of 
quantum numbers, and the operator C has the form 

C = T + V + 2Vh - ^i: (8) 
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where we again restrict our attention to stationary condensates. (The generahzation to the case of nonzero superfluid 
velocity is not difhcult In these Bogoliubov equations, the minus sign in the coupling terms is conventionally 

chosen to ensure that the ratio of the two amplitudes is positive for a uniform Bose gas, where plane waves are the 
appropriate eigenfunctions: u^e*'^ "" and Vfee**^ "". These coupled Schrodinger equations for the amplitudes Uj and Vj 
are analogous to the multicomponent Dirac equation. (The eigenvalues of the Bogoliubov equations come in ± pairs; 
the eigenfunctions of these pairs are related by a simple symmetry relation [^.) 

For a localized trapped condensate in an unbounded confining potential, the condensate density no and the two- 
particle Hartree potential Vr both vanish at infinity, and the amplitudes obey the usual quantum-mechanical bound- 
state boundary condition that Uj and Vj vanish for r — > oo. Furthermore, for positive eigenvalues, the eigenfunctions 
can be chosen to satisfy the orthonormality condition |9(] 



d r {u*Uk - v*Vk) = Sjk- (9) 

The (low-temperature) total noncondensate number N'{T) is obtained by summing over all eigenstates 

N'{T)=Y! N;{T), (10) 
j 

where the primed sum omits the lowest eigenstate (which simply describes the condensate itself). Here Nj{T) is the 
spatial integral of the corresponding temperature-dependent noncondensate density 

n;(r) ^ \v,ir)\' + {\u,ir)f + \v,ir)f){e^^^ - (11) 

The second term vanishes as T — > 0, leaving the integral of first term as the mean occupation number at zero 
temperature: 

7V'(0)= / dV|«,.(r)p. (12) 



At nonzero temperature, the condition No{T) = N — N'{T) determines the total condensate number, but the present 
work emphasizes the zero-temperature limit. 

It is helpful to recall briefly the special case of a uniform bulk condensate with constant density no [p|p^ , where 
the confining potential V is absent and the Hartree interaction energy and condensate chemical potential are equal: 
Vh = H = gnQ. The eigenfunctions are plane waves cx e*'' '", and the energy eigenvalues have the familiar Bogoliubov 



form 



hsk, for fc^ <C 1, 

h^P/2m, for fcC > 1, 



Ek = \lmvH + n^{ ^ ' (13) 



where = Ti^k'^ /2m is the kinetic energy and s — \J A'KaTi^nQ/m? — hj \f2m^ is the speed of compressional sound. 

The corresponding "coherence factors" and determine the mixing of the two components, and obey the 
bosonic normalization condition u^, — v\ = \ for each k. At zero temperature, the noncondensate occupation of the 
kth plane-wave mode is simply 

2 \(Tk^VH ^\ _ f v/TV87^« (2\/2A:e)-\ for fc^ « 1, 



' 2V ) \vllATl^\{kir\ forfc^»l. 

The long- wavelength singularity in A^^(O) cx is integrable, and the short-wavelength behavior ensures that the 
total noncondensate density n' — (27r)~^ J (fikv^ — |no\/noa^/7r is not only finite but also small relative to the 
condensate density. At low temperature, the corresponding additional thermal occupation for the kth state has the 
form AN'.iT) ^ N^T) ~ N',{0) (x T/k\ 



III. EQUIVALENCE WITH HYDRODYNAMIC FORMALISM 

We now return to the general inhomogeneous case. To proceed, it is useful to note that the GP equation (||) for the 
condensate wave function can be rewritten with Eq. (||) as >C^' = Vf/^', which suggests the following transformation 
of the Bogoliubov amplitudes 
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In particular, it is straightforward to verify that 



and 



Hi 

^/iVn 



{f + VH)U„ 



where T is a differential operator defined by 



Tf 



2m|*|2 



V- (|vl/pv/) = 



2m uq 



V- (noV/). 



To simplify (Q) further, define 



and 



Then the Bogoliubov equations can be rewritten exactly as 



(15) 
(16) 

(17) 

(18) 

(19a) 
(19b) 

(20) 



{T + 2Vh)G, ^EjF,. 
These two equations can be combined to give a single equation for Gj : 

{f'^ + 2fVH)Gj = E^Gj. 

The corresponding Fj follows from Eq. ( [l9b| ). Comparison with Eq. (H) shows that the normalization for the jth 
eigenstate is simply 

^ = ^Jd'r {F;G, + F,G]) -jfjd'r ^{F*G,). (21) 

The zero-temperature occupation of the jih excited state is then [compare Eq. (p])] 

N;^iO)^^Jd\\m'\F,-G,\\ (22) 

To make contact with the hydrodynamic description of these same normal modes, recall that the second-quantized 
operators for the fluctuations in the density p' and velocity potential are simply linear combinations of the field 
operators (j> and (p^ ||l^ . It follows that the corresponding normal- mode amplitudes pj and are also linearly related 
to the Bogoliubov amplitudes Uj and Vj. A straightforward comparison yields 



''-Wo 



Conversely, if the hydrodynamic amplitudes are known, the corresponding Bogoliubov wave functions become 



2im $, 



Pj im^ <I>j 
2*" ^ h 

Pj im'^ <I>j 
2* h 



(23a) 
(23b) 

(24a) 
(24b) 



A straightforward combination of Eqs. (|T^) and (|2^) immediately reproduces the known hydrodynamic equations for 
the normal modes of a stationary condensate ||Tl|-|13(|. For example, the product VnGj is just the density-fluctuation 
normal-mode eigenfunction pj itself (apart from a constant factor) and Eq. ( pO[ ) becomes 



V • (710 Vp,) 



4™^ 



V • <^ noV 



no 



_V • ( noV ^ 



no 



(25) 



where ujj — Ej/h is the normal- mode eigenfrequency. 
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IV. LOW-LYING NORMAL-MODE HYDRODYNAMIC AMPLITUDES IN TF LIMIT 



It is convenient to rewrite the Bogoliubov equations in terms of suitably rescaled variables. We let x = r/i?o be 
the dimensionless position vector, where Rq is the characteristic condensate radius, and introduce the dimensionless 
condensate wave function 

o3\ 1/2 

4.^0^) % (26) 

where the quantity in parentheses is roughly the volume per condensate particle. The rescaled condensate wave 
function x satisfies the dimensionless radial normalization condition J^x'^dx\x\'^ = ?70) where the dimensionless 
parameter fjo in Eq. ( p6| ) is defined by 

In the Thomas-Fermi limit, fjo becomes independent of the size of the condensate. To complete our scaling of variables, 
we express all energies in units of hajQ. 
The three relevant operators then become 

^ Vx-(lxpVx), (28a) 



2i?2|x 

y = ii?2x^ (28b) 

Vh^R^\x\'. (28c) 

These expressions clearly exhibit the dependence on the large parameter and show that the trap and Hartree 
energies are comparable to one another and both much larger than the kinetic energy. In the rescaled variables, the 
product TVh is independent of R, and the basic eigenvalue Eq. ( |2^ ) for the hydrodynamic normal-mode amplitude 
Pj becomes 



Vx-(|xl'V,p,) + ieVx-||xpV, 



\xr V Ixr 



(29) 



where Ej is the dimensionless energy (or frequency) of the jth normal mode, and e = i?^"* is the appropriate small 
expansion parameter. 

For a spherical trap, Baym and Pethick |^ showed that 770 ~ in the TF limit, and the rescaled condensate wave 
function is simply 

\x\'^xl^la~^')0{l-x), (30) 



where x is the scaled radial variable. For a spherical trap, the eigenfunctions of Eq. (^9|) can be written as a product 
of (real) radial functions Pni{x) and a spherical harmonic Ylm{0^(j))^ where j — {nlm), and n is the radial quantum 
number. 

Stringari has solved Eq. ( p9[ ) in the TF limit {i.e., to zeroth order in e) and has shown that the eigenfunctions 
and eigenvalues are independent of R, as discussed in Sec. I. In particular, the TF energy eigenvalue is given by 

Eli^l + n{2n + 2/ + 3) = a- i + 2n{n + a + 1). (31) 

where a = I -\- ^ is half an odd integer. The corresponding radial functions Pni{x) have the form x^ Pni{x'^), where 
Pniix^) are nth-order polynomials in x^. It is not difficult to verify that these polynomials satisfy the hypergeometric 
equation |23| with the explicit (unnormalized) form 

PM = F{-n, n + a + l;a + l;u)= u''^ — [ u"+'^ (1 - 7.)" ] ; (32) 

L [n + a + I) au"- 
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the first two such polynomials are Pqi = 1 and Pu = 1 — {a + 2)u/ (a + 1) = 1 — (2Z + 5)m/(2Z + 3). 

The polynomials Pni are a special class of Jacobi polynomials, with Pni{u) cx Prl"''^'' (1 — 2u). For any Z, they form an 
orthonormal set on the interval < m < 1 with weight m". An n-fold integration by parts with the explicit differential 
expression in Eq. (p2|) readily yields the radial normalization integral 



dx \pni{x)\ 



duu'' {PM\ 



r(a + l)n! 

r(n + 



1 



1 



2(2n + a + 1) ' 



(33) 



which also follows directly from the properties of the Jacobi polynomials . 

The transformation introduced in Eq. (^) explicitly eliminates the chemical potential from the Bogoliubov equations 
for the amplitudes Uj and Vj. As shown above, it leads to the hydrodynamic Eq. ( p9| ) that involves only the condensate 
density. In principle, the Bogoliubov description is wholly equivalent to this hydrodynamic description, but an explicit 
verification involves higher-order terms in the condensate wave function ||25|] . 

To understand the situation in more detail, recall that, in the grand canonical ensemble at zero temperature, the 
chemical potential /i determines both the total number of particles N and the number A^o < A^ in the condensate 
p^ . In the present approximation of retaining only the condensate contribution (i. e. N No), the Gross-Pitaevskii 
equation wholly characterizes the resulting functional dependence A^o(m)- For large dimensionless /i, this relation is 
just the familiar TF result plus small corrections |^^. It is straightforward to verify that direct substitution of the 
TF chemical potential and the TF wave function xo from Eq. (30) into the Bogoliubov equations yields an incorrect 
eigenvalue spectrum. A more careful treatment that includes the leading correction xi of order indeed reproduces 
the hydrodynamic eigenvalues and eigenfunctions. 



V. NONCONDENSATE OCCUPATION AT LOW TEMPERATURE 



Stringari's theoretical prediction of the frequencies of the lowest hydrodynamic normal modes rapidly received 
experimental confirmation |^,^. The present work extends his results to determine the occupation number of the 
same low-lying normal modes. This analysis makes essential use of the quantum-mechanical Bogoliubov amplitudes 
Unim and Vnim for quasiparticle creation and annihilation operators associated with these particular single-particle 
states and exploits the equivalence between the hydrodynamic and Bogoliubov descriptions. 

Equation (|2^) expresses the radial amplitude Gni directly in terms of pni', apart from a constant factor, we have 
Gni{x) CX p„i(a:)/|xo(a;)P, where |xo(2^)P = ^(1 — x^) is the parabolic TF condensate density profile. In addition, 

2i?2i 



Eqs. ( |l9bD and (g8^) show that Fni ~ 2i?^|xorG„; apart 
Together, these expressions suggest the following normalization 



from corrections that become small as R 



Fnl{x) 



2R Cnl 
Enl 



Pnl{x) 



2RCnl I f 2\ 



■x'Pnlix'), 



(34a) 



Gnl{x) 



Cnl Pnl{x) 



R \xo{x)\^ 



Cnl x'Pnljx^) 

R \xo{xW ' 



(34b) 



where Cni is a normalization constant determined from Eq. (21). A combination of the previous results leads to the 
explicit expression 



^2 _ 



27rTyo Enl 

TO 

^nl 



(35) 



which completely determines the radial amplitudes Fni and Gni in the TF limit. Note that while Fni is of order R, 
Cnl is of order 1 /R, so that Uni = \ {Fni + Cni ) and Vni = ^ {Fni — Cni ) are both large for the low-lying states of a 
large condensate. [A similar behavior occurs at long wavelengths {k^ ^ 1) for a uniform dilute Bose gas, as seen from 
the bosonic normalization condition Uk ~ ~ ^ Fq. (|l^) .] 

The zero-temperature occupation of the low-lying normal modes with quantum numbers nlm follows directly from 
Eq. (|2^). The leading term from \Fni\'^ is of order i?^, and the cross term between Fni and Cni just reproduces the 
normalization integral from Eq. (p3) . It is convenient to introduce the general class of integrals 



F 

^nl 



''dx{l 



[P7U{X)]' 



duu'^il-uy [Pnl{u)y 



(36) 
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For j = this quantity is simply the normahzation integral considered previously; for j 
repeated integration by parts, leading to the ratio 



1 , it can be evaluated with 



2n{n + a + 1) + a 
{2n + a){2n + a + 2) " (2n + a)(2n + a + 2) ' 



(37) 



A combination of these results yields an explicit expression for the zero-temperature occupation of the low-lying 
excited states of a large isotropic condensate in an isotropic harmonic trap: 



AEni (2n + a)(2n + a + 2) 



1 

2 ' 



(38) 



where the omitted term comes from the radial integral of |G„ip and is of order . 

Since i?^ ^ 1 in the TF limit, the noncondensate occupation N^i{0) is large for the low- lying states with small n 
and I, whose dimensionless energy is of order unity. For example, for the lowest dipole mode with n — 0, I — 1, we 
have Eoi = 1 and A^oi(O) ~ T4^^ ~ \- "^^^^ behavior is evidently very similar to that for a uniform dilute Bose gas, 
where Eq. (|l|) shows that iV((0) = {2^/2k^y^ > 1 for long wavelengths such that < I. Indeed, Eq. (^ also 
makes clear that the long-wavelength approximation fails when the kinetic energy Tk of the plane-wave state becomes 
comparable to the Hartree energy Vh, and a similar behavior is expected in the present case of a trapped condensate. 

Specifically, the amplitude F^i was obtained from Gni by neglecting T relative to 2Vh in Eq. ( [19b| ). This approxi- 
mation holds for radial states with sufficiently few nodes, but it necessarily fails for highly excited states whose large 
kinetic energy reflects the bending energy associated with rapid oscillations (and hence many nodes) in the wave 
function [Tl| ]. By analogy with the corresponding situation for a uniform Bose gas, it is natural to conjecture that 
the explicit expression in Eq. (38) is valid only for low- lying modes with iV^;(0) ^ 1. Verification of this conjecture 
would require a detailed study of the highly excited modes of the large condensate; it involves corrections to the 
TF condensate wave function associated with the boundary layer ||2l[ and a WKB (phase-integral) description of the 
rapid oscillations inherent in the short- wavelength limit. This difficult analysis remains for future investigation. 

The total number TV' of noncondensed particles is the sum of iVj over all eigenstates of the Bogoliubov equations, 
omitting the lowest solution with zero energy (which describes the condensate itself). In the present case of a spherical 
condensate at zero temperature, we have 



N'{o) = Y^ (2; + iX,(o), 



(39) 



nl 



where the factor 21 -\- 1 represents the degeneracy associated with the sum over azimuthal quantum numbers m and 
the prime on the sum indicates that the term n = ^ = is omitted. In the TF limit, the zero-temperature occupation 
of the low- lying modes is given by Eq. (^). As argued above, the sum must be cut off when N^i{0) « 1. (Apart from 
a numerical factor of order unity, an analogous cutoff gives the correct total noncondensate fraction for a uniform 
condensate at zero temperature, N'/Nq y/noo^). Stringari has suggested that the TF expression for Eni 
in Eq. (|3^) holds for E„i Mi which is ^i?^ in the TF limit. It is not hard to see that this criterion provides a 
qualitatively similar cutoff. 

For large R, the N^^ vary slowly with n and and the double sum Eq. (39) can be approximated by an integral 
over continuous variables n and a = I + ^, with 



N\0) = y(2l + l)KiiO) 



71 1 



dn / da 2a 



r:' 



AEni {2n + a){2n + a + 2) 



(40) 



where this double integral runs over the region N'^^i 
variables 



> 1. To clarify its structure, it is convenient to introduce new 



s = n-|-/-|-l = n-|-Q;-|-i 



and 



(41) 



or, equivalently 



and 



(42) 



this transformation has unit Jacobian, and the allowed region in the st plane is the first octant < i < s, apart from 
a small region around the origin. In these new variables, the TF energy eigenvalue in Eq. ( pl| ) becomes E"^ = 2st— 1, 
and the corresponding noncondensate occupation number is 



8 



KM 



4 



2st 



^2st - 1 (s + t)2 - 1 



(43) 



To isolate the dominant contribution to the integral in the st plane, it is helpful to introduce plane-polar coordinates 
(C, (p), with s = Ccos0 and t — (siiKp; for large the leading behavior is N' ~ i?^/C, apart from angular factors. It 
is not difficult to see that the angular integral over < </> < 7r/4 converges, and the radial integral must be cut off at 
Cmax ^ R^- In this way, the total noncondensate number at zero temperature becomes 



N'{0) 



c 



(44) 



As shown by Baym and Pethick |8 
proportional to R^. So from Eq. (p 
Thomas-Fermi limit is 



the number of condensed particles in the TF limit is given by A^o ~ j^doR^ /a, 
h, the ratio of uncondensed to condensed particles at zero temperature in the 



N'jO) 
No 



do 



(45) 



apart from a numerical constant of order unity. As expected, this is comparable to the uncondensed fraction in a 
homogeneous Bose gas of density Nq/Rq. Since the validity of the Bogoliubov approximation depends on the condition 
N' /Nq <C 1, we find that the present description holds only for R <C do/a. When combined with Baym and Pethick's 
result, this condition indicates that the condensate is dilute only for No ^ R^{do/a) ^ (do/a)^, in agreement with 
the condition found in Sec. II. 

It is evident from Eqs. (|) and (|ll|) that the zero-temperature eigenvalues and eigenfunctions of the Bogoliubov 
equations also determine the low-temperature thermal depletion of the condensate through the relation 



1 + 2A^;K0) 
exp{phuJoEni) - 1 ■ 



(46) 



where Eni = \/l + n{2n -I- 2^ -I- 3) is the dimensionless energy eigenvalue, here taken from Eq. (^l]). It is convenient to 
define a characteristic temperature O = haJo/kg for thermal excitation of the first excited state; note that Q is much 

smaller than the ideal-gas transition temperature Tq^'' ~ 9A^^/"^. Furthermore, the actual transition temperature To 

for a dilute trapped Bose gas is only slightly less than the ideal-gas transition temperature Tq^\ so that 9/To « N^^/'^. 
As a result, the low-temperature thermal occupation of each low-lying mode increases linearly with temperature 



T 



eEni 



(47) 



for OEni < r< To. 



VI. DYNAMIC STRUCTURE FACTOR S(q,tj) 

Consider an external probe that scatters with momentum transfer hq and energy transfer huj to a target. If the 
probe couples weakly to the number density of the target (here, the trapped Bose condensed system), the differential 
cross section d^a/dfldcu is proportional to the dynamic structure factor [ p6| , p7[ 

1 ^ .-l3E,,,r,~nu\\2:f,„^^-^^' 



5(q, -) = ^ 1^ e-^-' \{f\p'^m^6lu;- 1 , (48) 

fi 

where i and / refer to exact states of the interacting target with energies Ei and Ef, Z = J^i exp(— /JE".;) is the target 
partition function, and — J d^r e^'^'^ /5'(r) is the "creation operator" for a density fluctuation with wave number q. 
In the Bogoliubov approximation, the density-fluctuation operator p'(r) is proportional to a linear combination of 



quasiparticle operators |13] a, and a\ 



the field operators (/)(r) and 0t(r). As a result, Eq. (23a) immediately yields the corresponding expansion in bosonic 



p'{v) = Y^[p,{v)a,+p*{v)a% (49) 
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where pj is essentially a linear combination of the Bogoliubov amplitudes, and the primed sum runs over all the 
excited states of the condensate. The evaluation of the dynamic structure factor is straightforward, giving the explicit 
result 

^(q'^) = ^ E' [(1 + Z^) IP;(q)l''5(^ - E,lh) + /, \PM)?5{u: + E,/h)\ , (50) 







where fj = [exp(/3i?j) — l] ^ is the thermal Bose-Einstein function and /5j(q) is the spatial Fourier transform of 

. . . ■ vr^ 

For a spherical trap, the amplitudes are given in Eq. ( p4b| ). Introducing the dimensionless variable Q = qRn, we 

readily obtain the following dimensionless dynamic structure factor (scaled with the oscillator frequency luq) 

S{Q, io) « b„KQ)l' [(1 + fnl) Sioj - E,u) + fnl Sioj + E,u)] , (51) 

nl ""-^ 

where Eni is the dimensionless energy (or frequency) from Eq. (^) and 

Pni{Q)= I dxx'+^P^i{x^)ji{Qx). (52) 



Equation (|5T| ) makes it easy to verify that this approximate dynamic structure factor obeys the detailed-balance 
condition gej S{-Q,~lj) = e-l^'''^S{Q,uj). 

The integral in Eq. (m) can be evaluated with the standard expression for the spherical Bessel function ji as an 
^fold derivative of jo [ |28[ 

^$a.,_2y^,„M, (53) 

where w = cc^, along with the explicit formula for P„i(x^) in Eq. (^2|): 

/n^ r(a + l) 2'-H-l)' f\ r „+„.^ . , 

^"'(^)= r(n + a + l) ^^ L (1-..) ]^jo(gx). (54) 



An 7i-fold integration by parts then yields 

r(a- 

r(n + a + l) 2«+i 7o (Qa;)"+' 

which can be integrated term- by-term after expanding the spherical Bessel function as a power series in uQ^. The 
resulting series can be re-summcd to yield the remarkably simple expression 

r(Q; + l)n! j2n+;+i(Q) 

Pnl[Q) = 7=^— — TT 7^ • 56 

1 (n -I- a + 1) Q 
It is convenient to rewrite the dynamic structure factor as 

SiQ,u) = Y^'SnliQ) [(1 + fnl) 6(10 - Enl) + fnl S{iO + E^i)] , (57) 
nl 

where /„; = [exp(e£;„i/r) - 1]"\ and 

SnliQ)^ ~^Q2Ji2 [j2n+l+liQ)\ ■ (58) 

For Q — + 0, the leading term arises from the dipole-sloshing mode n = 0, I = 1, with ^oi « /2R^ + 0{Q^); the 
next contributions (of order Q^) arise from the terms with n = 1, I = and n = 0, I = 2. In addition, the frequency 
integral is simply the static structure factor 

cLoSiQ,u;) = 5(Q) 5„i(g) coth j , (59) 
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with the long-wavelength limit 



5(Q)«^COth(^|;)+O(04). 



(60) 



Finally, the first moment is the /-sum rule | |2^ 




duju) S{Q,u;) = Sni{Q)E, 



( 




in conventional units) . 



(61) 



2™ 



The dipole-sloshing mode exhausts this sum rule at long wavelengths, but it otherwise implies a rather intricate 
identity involving sums of squares of spherical Bessel functions; it is easy to verify this relation through terms of order 
Q^, but we have not sought an independent derivation. 



This work has shown how the known hydrodynamic amplitudes pJ| | for a dilute condensed Bose gas in a spher- 
ical harmonic trap can provide the corresponding Bogoliubov spatial amplitudes for the quantum-mechanical field 
operators, including their absolute normalization. These normalized amplitudes in turn determine several important 
physical quantities, such as the excitation of each normal mode, both in the ground state and at low temperatures 
(T <^ To), and the dynamic structure factor that describes the inelastic scattering of a weakly interacting probe 
that couples to the density fluctuations of the Bose condensate. The hydrodynamic amplitudes are accurate only 
for low-lying modes, for their derivation neglects the quantum-mechanical kinetic energy that becomes increasingly 
important at short wavelengths. Thus our estimate of the total condensate depletion involves a cut-off that explicitly 
omits all the high-lying normal modes. 

As noted by Singh and Rokhsar |^ and Stringari |l^] , the exact eigenfrequencies of the lowest dipole mode ^ = 1 for 
an anisotropic harmonic trap coincide with the bare oscillator frequencies; numerical work by several groups ]29|-|3^] 
has confirmed this conclusion for moderate values of Nq. The Appendix contains an analytical proof in the Bogoliubov 
approximation; it provides an explicit construction of the lowest dipole states for any solution of the GP equation, 
including those containing one or more vortices. 



We are grateful to S. Bahcall, M. Cole, A. Grifhn, and S. Stringari for valuable comments and suggestions. This 
work was supported in part by the National Science Foundation, under Grants No. DMR 94-21888 and DMR 91-57414. 



Rece ntly, t he collective excitation spectrum of harmonically confined Bose gases has been calculated by several 
groups ||ll| , |29| -|33|| . It has been noted ||ll|,^ that these spectra contain a trio of exact collective modes that correspond 
to the simple-harmonic oscillation of the center of mass of the condensate. If the interacting condensate is displaced 
without deformation, the interparticle interactions are unchanged. Thus, the system experiences a restoring force 
that is simply linear in the displacement, with spring constant equal to that of the bare trap. There is one such mode 
for each of the three Cartesian directions. 

In this Appendix, we construct the exact many-body raising operator that creates these "dipole" or "sloshing" 
modes and demonstrate the independence of their frequencies on the interparticle interactions. Since the collective 
modes of the dilute Bose gas are expected to be well-described by the Bogoliubov approximation, it is interesting 
to confirm that the dipole modes are in fact unrenormalized in this approximation. We therefore construct the 
exact solutions of the Bogoliubov equations that correspond to the dipole modes. These may be useful points of 
comparison with numerical calculations of the Bogoliubov spectrum. Finally, we examine the form of these modes in 
the Thomas-Fermi approximation which is relevant for large condensates. 



VII. DISCUSSION 
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APPENDIX: DIPOLE MODE 
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Definitions. 



Consider bosons of mass m confined to a three-dimensional, harmonic potential with spring constants muj'^, where 
a = X, y, and z labels the three cartesian directions. The corresponding bare trap frequencies are uja- The (one-body) 
Hamiltonian of a single particle in the trap is then 



r,2 



^o^E(£ + l--Mj, (Ala) 
or, equivalently, 

Ho=Y, ?i^a(4aa + 5), (Alb) 

a 

where we introduce the usual raising and lowering operators 



\/2 ( c?Q " dxa ) ' ^ ^ 



and (as above) da = \ff\JrrvjJa. The operators and a]j obey Bose commutation relations, and they act on one-body 
states. They satisfy the familiar relation 

[Ho,aJ,] = TiWaaJj (A3) 

which shows that is a raising operator of the noninteracting system. 

Let V(v — r') be the two-body interaction potential between the particles. The many-body Hamiltonian is then 

N N 

H = EHo. + -Enr.-r,), (A4) 
j=i i^j 

where the index i and j labels individual particles, and N is the total particle number. For bosons, we consider only 
states that are symmetric under interchange of particle labels. 



Dipole-mode creation operator 

We claim that the (symmetric) operator 



1=1 

is a raising operator for the many-body Hamiltonian 7i. That is, 

[U^AlX^hujaAl. (A6) 

Thus if |G) is the exact many-body ground state of the interacting system, then AjjjG) is an exact excited state with 
excitation energy husa- Repeatedly applying to any exact eigenstate of 7i creates an equally spaced ladder of exact 
excited states. [Note that (A„)^ creates two quanta of the same elementary excitation and does not yield a distinct 
new elementary excitation itself.] 



The proof is straightforward. The commutator Eq. (A6) is 



[H,At] = + lEt^'^'^iJ- (A7) 

ii ijk 

The first set of commutators is easily evaluated using Eq. ( |A3| ) and gives ?iWa X^i ^L— ^^"^q- The second set of 
commutators vanishes identically, since the terms cancel in pairs. 

This operator A]^ makes good physical sense. Beginning from any exact many-body eigenstate, we raise each particle 
in turn by a single quantum, and then superimpose the resulting states to produce a symmetric wave function. The 
proof is similar to the demonstration that the cyclotron frequency of an interacting, translationally invariant system 
is unrenormalized, i.e., Kohn's theorem [p4l. 
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Dipole mode in the Bogoliubov approximation. 



Is the same dipole mode also present as an exact solution of the Bogoliubov equations for harmonically trapped 
bosons? The answer to this question is not completely obvious, since the Bogoliubov approximation assumes the 
existence of a condensate, which provides a provides a preferred reference frame (usually, but not necessarily, at rest). 
We will show that the spectrum of the dilute Bose gas in the Bogoliubov approximation possesses exact excited states 
that are simply related to the dipole modes discussed above. 

Let ^'(r) be the condensate wave function of the interacting system, which is an exact solution of the Gross-Pitaevskii 
equation [Eq. (||)]: 



(7io + g|*n* = M*, 



(A8) 



where g — Anah^ /m characterizes the strength of the interparticle potential and fi is the chemical potential. We claim 
that the excited state 



Ua{r) 



4*(r) 
aa**(r) 



(A9) 



is then an exact solution of the Bogoliubov equations with excitation frequency equal to the bare trap frequency uja ■ 
There is one such mode for each coordinate direction. 

Note that for the noninteracting case {g = 0), annihilates the condensate wave function which is then simply 
the ground state ^'o of the harmonic potential. Equation (A£) then reduces to 



Ua(r) 
Wa(r) 



4*o(r) 




(noninteracting) . 



(AlO) 



As expected, exact quasiparticle states are created by adding a particle to the first excited state of the trap. 

It is straightforward to prove that Eq. ( [A9| ) is an exact solution of the Bogoliubov eq uati ons for any coupling 
strength g. The proof proceeds by the explicit demonstration that the Mq, Va of Eq. (A9) satisfy the coupled 
Bogoliubov equations (0). We use the elementary properties of Oq and aj,, and the fact that satisfies the Gross- 
Pitaevskii equation [Eq. ( |A8[ )]. The argument explicitly relies on the harmonic character of the trap V{v), with aj, 
and Cq, as raising and lowering operators for the single-particle problem; it cannot be extended to problems without 
a simple ladder of noninteracting levels. 

We now confirm that the first of the two coupled Bogoliubov equations Eq. (^ is satisfied by our proposed 
solution EcL (A9); the second equation Eq. (^) is simply the complex conjugate of the first. Substituting Eq. ( [A^ ) 
into Eq. (|7a|), we need to show that 



[Ho + .g|*|' - /i) 4* + .g|*|'al* - g^^a^^* I huj^al^^. 
Multiplying the Gross-Pitaevskii equation (^ by aj, yields 

al(7^o+g|*P-M)* = 0. 



Subtracting Eq. ( A12 ) from Eq. ( All ) then reduces our problem to showing that 

[Ho, aU^ + .g[ 4]* -I- gl^l^flt ^ - 5*2^^^* I hLOo.ai^. 



(All) 
(A12) 

(A13) 



From Eq. (A3), the first commutator is simply [Ho, aj,] = Tiuja^a^ which cancels the right hand side of Eq. ( Al3| ). 
Finally, our task is reduced to showing that 



(A14) 



this result is easy to confirm with the explicit form of a\ and Cq, in Eqs. ( pa| ) and ( |A2b| ). 

Note that this construction is entirely general. It makes no assumption that 5* is real, and thus holds for any solution 
of the Gross-Pitaevskii equation (||), including those describing vortices [^. For any self-consistent condensate, the 
lowest dipole modes will have the bare oscillation frequencies. 
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Nature of the dipole mode. 



Excitation of the dipole mode leads to an oscillatory density-fluctuation amplitude that is one of the normal modes 
described in Eq. (p3|). For simplicity, we consider explicitly only the case of a condensate with real i.e., a stationary 
condensate. 

Equations ( p3| ) then become 

= V-2^d^^ocp^. (A15) 

This expression confirms that the quasiparticle mode created by the linear combination of aj^ and corresponds to 
a uniform displacement of the condensate in the ath direction. Similarly, the velocity potential associated with this 
condensate motion is [see Eq. (p3|)] 



_ ~iT^ Tg _ (A16) 
\/2m da 

the corresponding velocity v = V(5<I>a lies a long the displacement with constant amplitude and is ^tt out of phase 
relative to the density fluctuation Eq. ( A15) due to the factor of i. 



Thomas-Fermi approximation 



Finally, we examine the dipole mode in the Thomas- Fermi limit, where the Gross-Pitaevskii equation can be solved 
exactly. For simplicity we consider an isotropic trap, although the same results hold for anisotropic traps of the sort 
used in Refs. |^,^. The condensate wave function is then given by 



where J^o oc N^^^ is the size of the condensate. Using Eq. (24) and the operators from Eqs. (A2), we find 

~ 1 



,(r) = aj,1'(r) cx 



(A17) 



(Al8a) 



Va{r) = aQ*(r) cx 



(A18b) 



Although dp and (5$ remain finite everywhere, both both and Va diverge at the perimeter of the condensate. 
This pathology evidently arises from the extreme Thomas- Fermi limit, for it reflects the singular behavior of "^tf 
at the condensate boundary. Any finite interparticle interaction strength orenders the condensate wave function ^ 
differentiable everywhere, and the exact 4" vanishes smoothly for r — s- cx3 then Ua and Va also vanish smo othly . 
Despite the separate divergences, the approximate two-component Bogoliubov state in the TF limit in Eq. (AIS) 
remains normalizable, since we require only that Jg dr — |waP) = 1- 
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